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Abstract. Using mid- infrared and radio criteria, we select a sample of candi- 
date z ~ 2 obscured quasars. Optical spectroscopy confirms about half of these 
as type-2 quasars, and modelling the population suggests 50-80% of the quasars 
are obscured. We find some flat radio spectrum type-2 quasars, and tentative 
evidence for obscuration unrelated to the torus. Using a similar sample, we also 
find evidence for a significant fraction of Compton-thick quasars. 



Introduction 



The population of high-redshift obscured (type-2) quasars (active galactic nu- 
clei with Lboi^ 10 38 W) has proven an elusive one, despite great progress in 
deep surveys. In radio-selected samples, the fraction of unobscured (type-1) 
quasar s to the total num ber of quasars (the quasar fraction) is found to be 
0.5 ( Willott et al.|[2000l ). corresponding to a ratio of type-2 to type-1 ~1:1. 



At X-ray wavelengths, surveys have been successful in finding large numbers of 
type-2 AGN, but the exact fractio n of type-2 quas ars is still unclea r (for recent 
results, see e.g. IWang et al.ll2006l : iDwellv fc Pagell2006l : iGilli et al.1 [20061 k 



A possible explanation for the discrepancies is that some quasars are very 
heavily obscured (with absorbing columns, -/VhkI 10 28 m~ 2 ) and are therefore 
very difficult to detect in X-ray surveys. An alternative approach is to use the 
characteristic warm dust (T ~ 100-1000 K) signature of AGN to select obscured 
quasars based on their mid-infrared properties, and indeed lLacv et ail (|2006l ) 
confirmed the success of a 3.6 to 8 /im colour selection to find AGN (obscured 
and unobscured) at intermediate redshifts. 

Of particular cosmological relevance is the quasar fraction at z ~2, around 
the maximum in the quasar activity. We therefore devised a set of strict selection 
criteria to look for z ~2 type-2 quasars. This has also practical advantages, since 
above z = 1.7 Lya is visible with ground-based optical spectrographs. 



2. Selection 

Quasars are bright at mid-infrared wavelengths due to the presence of warm 
dust. At z ~2, observed 24 /im corresponds to rest-frame ~ 8 /xm, or dust 
at ~ 360 K. At these wavelengths, the obscuration due to dust is negligible, 
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unless the Ay is extreme. The extinction due to dust increases rapidly with 
decreasing wavelength, so that type-2 quasars that are bright at 24 /im should 
be significantly fainter at 3.6 /im (rest-frame ~1.2 jttm). In addition, if the 
quasar has a value of Ay^, 10, the 3.6-/mi flux density will be dominated by the 
host galaxy, instead of the quasar light. Since quasar host galaxies are typically 
massive ellipticals, one can assume a typical galaxy mass, and use the 3.6-fim 
flux density to obtain a crude photometric redshift. 

Therefore we demand objects to meet the following criteria: S24 ^ m > 
300 fiJy and S3. 6 < 45 /iJy Assuming a typical quasar SED, the 24-/mi 
criterion corresponds to slightly below the break in the quasar luminosity func- 
tion at z = 2, while the S^^fim flux density cut corresponds to z p h ot > 1-4. This 
allows some margin of error for our crude z p h ot cut. 

These two criteria should select quasars with Ay^ 5 and Zphot^ 1.4. How- 
ever, they are also likely to include a large number of ultra-luminous infrared 
galaxies (ULIRGs) at z ^1. We therefore included a radio criterion: 350 /iJy < 
Sia GHz < 2 mJy. Only ULIRGs with very high star-formation rates 100-350 
Mq yr _1 of massive stars) are expected to make this cut, which therefore min- 
imises the contamination in our sample. The trade-off is that the quasars we 
select are radio-intermediate (Li.4GHz~ 10 24 WHz _1 sr~ 1 ), so although their 
(rest-frame) 8-/im luminosity densities correspond to below the break in the 
luminosity function, their radio properties make them more rare than the radio- 
quiet quasars. 



3. The quasar fraction 

Applying this criteria, we found a number of candidates (21 in the 4 deg 2 of the 
Spitzer extragalactic First Look Survey, FLS), and optical spectroscopy con- 
firmed about half of these as narrow-line type-2 quasars at z > 1.7. One object 
at z = 1 had no high-ionisation lines, and could plausibly be a ULIRG contam- 
inant. The rest of the objects had co mpletely blank spectra, and the re is no 
contamination from type-1 quasars (see Martmez-Sansigre et al.| [2006a). 



To estimate the quasar fraction, we modelled the expected number of type- 
1 quasars meeting our 24-^m and 1.4 GHz criteria, and at z > 2.0, where our 
completeness is expected to be high. Assuming a flat prior probability distri- 
bution (PDF) for the quasar fraction (reflecting our ignorance), we calculated 
the posterior PDF of the quasar fraction, given the predicted number of type- 
1 quasars, and the number of type-2 quasars showing narrow- lines and with 
z spec > 2.0. This is shown as the solid curve in Figure [TJ The distribution peaks 
around a value of the quasar fraction, q, ~ 0.5 implying a similar number of 
narrow-line type-2 quasars as type-1 quasars. 

Some of the blank objects will be in the "spectroscopic desert" (1.2 < z < 
1.7), yet many are probably at z > 1.7 (and two cases have been spectro- 
scopically confirmed to be at z = 1.75 and 2.25 from Spitzer spectroscopy by 
IWeedman et al.l [20061 ). We therefore use the crude photometric estimates (from 
the 3.6 /im flux density) to estimate how many blank-spectrum objects lie at 
z > 2.0. When these objects are added to the narrow-line objects, one obtains 
the dashed-line distribution of Figure [H and the modal value for q is now closer 
to 0.25. The true quasar fraction is probably somewhere between q = 0.2 and 
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Figure 1. Normalised posterior probability distribution of the quasar frac- 
tion, given the data and the predicted number of z > 2.0 type-1 quasars 
matched in 24-/im and radio criteria. The solid line shows the distribution if 
only the narrow line type-2 quasars at z S pec> 2.0 are taken into acount, while 
the dashed line shows the distribution if the blank objects with z p hot> 2.0 are 
also included. 



q = 0.5, so that we can conclu de that at the peak of type-1 quasar activity, 
50-80% of quasars are obscured ( Martmez-Sansigre et al.l 120051 ). 

A possible explanation for the blank objects is that they are being obscured 
by dust on scales larger than the narrow-line region. The obscuration would 
therefore be independent of the orientation. If this scheme is true, then some 
of the blank-spectrum objects are expected to harbor quasars with the torus 
face-on: objects that would be type-Is if they were not obscured by dust on kpc 
scales. If the radio jet is close enough to our line of sight, such a quasar might 
exhibit a flat-spectrum. 

To test this, we used radio data at 1.4 and 4.9 GHz from the VLA and at 610 



MHz f rom the GMRT to measure the radio spectral indices ([Martmez-Sansigre et al. 



2006b). Most objects show steep spectra, as expected from radio sources, and 



some gigahertz-peaked sources are also found. However, four objects show flat 
radio spectra. One of these is the z = 1 object, and its radio luminosity 
(■^i.4GHz~ 10 23 WHz _1 sr -1 ) is consistent with a massive-star-formation rate 
of ~300 Mq yr . it is probably a ULIRG contaminant. The other three ob- 
jects are probably flat-spectrum type-2 quasars, and none show narrow-emission 
lines (one of them has z spec = 2.25 from Spitzer spectroscopy). 



4. Compton-thick quasars 

If this scheme is correct, as suggested by the evidence collected so far, then type- 
2 quasars outnumber type-1 quasars by ~ 2-3:1, and therefore a fraction of type-2 
quasars seems to have been missed out by hard X-ray surveys. The most likely 
explanation is therefore that some of these type-2 quasars are optically-thick 
to Compton-scattering, i.e. are "Compton-thick" (Nn> 10 28 m~ 2 ). Since no 
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X-ray data was available for the first sample, we selected a similar sample in the 
Subaru-XMM Deep Fiels (SXDF) where deep X-ray data was available. The area 
here is smaller (0.8 deg 2 ), so we needed to change our selection slightly to yield 
more quasars. We lowered the lower cut in radio flux density to Si.4GHz > 100 
//Jy and left all other criteria unchanged. As well as increasing the number 
of high-redshift type-2 quasars, this would increase the number of ULIRG con- 
taminants, and no spectroscopy was originally available to distinguish between 
objects. 

To separate the type-2 quasars from the ULIRGs, we fitted the spectral 
energy distributions (SEDs) between B-band and 24 jum, using a quasar and a 
ULIRG model, and used a bayesian criterion (the odds ratio) to select between 
models. We only kept the objects where the quasar model was strongly favoured 
over the ULIRG model, with reasonably good fits, and with Zp^Qi > 1.70, so that 
they could be tested observationally. This left us with 12 good type-2 candidates, 
only two of which were detected in the deep X-ray catalogue of the SXDF. From 
the photometric redshifts and the X-ray data, the two objects have estimated 
column densities of logio(A r H/m 2 ) = 27.25 and 27.50, so they are both heavily 
obscured but Compton-thin. 

The other 10 objects are undetected in X-rays, up to an approximate flux 
limit in the 2-12 keV band S2-i2keV~3x 10 -18 W m -2 . They have similar es- 
timated values of z p hota n d bolometric luminosities to the two X-ray detections, 
so if they are indeed quasars, they are probably more heavily absorbed, with 
A^h^ 10 28 m -2 . We therefore find evidence that many of these objects are 
Compton-thick (Martmez-Sansigre et al. 2006c). 

Acknowledgments. We thank our collaborators for their important con- 
tributions to this work. 
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